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The recent Cosmic Microwave Background (CMB) measurements indicate that there is power 
deficiency of the CMB anisotropics at large scales compared with the ACDM model. We have 
investigated the possibility of explaining such effects by a class of primordial power spectra which 
have infrared cutoffs close to the horizon scale. The primordial power spectrum recovered by direct 
deconvolution of the observed CMB angular spectrum indicates that the data prefers a sharp infrared 
cutoff with a localized excess (bump) just above the cutoff. We have been motivated to assess 
plausible extensions of simplest inflationary scenarios which readily accommodate similar form of 
infrared cutoff. We carry out a complete Bayesian analysis of the parameter space using Markov 
Chain Monte Carlo technique with such a class of primordial power spectra. We show that primordial 
power spectrum that have features such as an infrared cutoff followed by a subsequent excess in 
power give better fit to the observed data compared to a nearly scale-invariant power law or power 
spectrum with just a monotonic infrared cutoff. However, there is substantial room for improvement 
in the match to data and calls for exploration of other mechanisms that may lead to infrared cutoff 
even closer to that recovered by direct deconvolution approach. 

PACS numbers: 98.80.Cq, 98.80.-k 



I. INTRODUCTION 



The Cosmic Microwave Background (CMB) is the most copious of all the isotropic backgrounds observed 
from radio to gamma rays. It has become the most important observational tool both for determining the 
global cosmological parameters and to understand the physics of the early universe. The concordant model 
of cosmology that has emerged from recent observations is the inflationary ACDM model. This describes the 
universe as spatially flat with primordial fluctuations which are adiabatic and Gaussian. The ACDM model 
has been quite successful in explaining many aspects of the universe conforming most with the observational 
data [1-4]. The primordial power spectrum formed during the inflationary epoch, in the simplest case can 
be described as a single inflaton field with a nearly scale-invariant Harrison-Zel'dovich spectrum expressed 
as V(k) oc A:™ 3-1 , 'fc' being the comoving wavenumber. However, inflationary scenarios readily accommodate 
significant departures from scale invariance [5-11]. Distortions in the shape of the initial power spectrum 
imprints characteristic scales on the form of V{k). The amount of such distortions manifested on the angular 
power spectrum and the scales at which they occur are strictly related to the nature of the primordial 
fluctuations and to the matter content of the universe. Careful analysis of all features of the CMB power 
spectrum put constraints on the entire parameter space. Recently a lot of attention has been devoted 
to the power anomalies at low multipoles [12] compared to that obtained using the ACDM model of the 
universe. This was first reported by the Cosmic Background Explorer (COBE) [13] satellite and now further 
confirmed by Wilkinson Microwave Anisotropy Probe (WMAP) observations [2, 14, 15]. In the recent past 
many speculations have been made to provide a satisfactory explanation to this discrepancy but they have 
been successful only to a certain extent. The WMAP team have crudely fit the power at low scales by putting 
a running spectral index to the nearly scale-invariant power law [14]. This aspect has also been projected 
as a probable new feature in the CMB anisotropies such as related to curvature scale [16]. This has been 
related to the topology of the universe as well, such as a dodecahedral universe [17] or possibly a violation 
of statistical isotropy [18]. A low quadrupole moment has been suggested to be originating also from dark 
energy, quintessence models with isocurvature perturbations [19, 20] via correlations with adiabatic density 
fluctuations which lower the adiabatic Integrated Sachs- Wolfe contributions to the quadrupole. Even model- 
independent searches have been made to look for features in the CMB power spectrum [21-25]. The studies 
where the power spectrum has been reconstructed in a set of bins in k-space have the obvious disadvantage 
that the reconstruction depends on the width and the placement of the bins. This is also true for methods 
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such as linear interpolation, wavelet, top-hat bins which can contribute to smoothening out or omission of 
features. Despite the progress the issue of low scale power still remains a subject of debate. The most 
crucial point perhaps is the limited information available at large angular scales due to cosmic variance [26] . 
Although large cosmic variance at the quadrupole does make the power discrepancy less significant, it is 
important to bear in mind that cosmic variance is proportional to the 'true' angular power spectrum of the 
underlying model. The cosmic variance is small for a model that predicts the observed low quadrupole. In 
principle, the case for an infrared cutoff is best addressed by bayesian evidence study. However, different 
bayesian evidence studies for the models with infrared cutoff has yielded different conclusions regarding the 
evidence for an infrared cutoff in the primordial spectrum [28, 29] . This is possibly due to the dependence of 
bayesian evidence analysis on the different class of models explored, details of parameterization of models and 
the datasets used. Hence these results may not readily translate into a general conclusion on the evidence 
for 'infrared cutoff' though [28] find that models with cutoff in power generally produce higher evidence 
compared to scale invariant or a single spectral index model. 

We revisit the possibility of explaining this observed suppression of power at low multipoles using known 
mechanisms of modifying the shape of the primordial power spectrum from inflation to introduce an infrared 
cutoff [30-33]. The limited success of model spectra with monotonic infrared cutoff studied is apparent 
from the non-monotonic infrared cutoff indicated by direct deconvolution of the primordial spectrum from 
the CMB angular power spectrum using WMAP data [24, 34]. The direct approach suggests that the data 
favours a sharp infrared cutoff followed by a localized excess ('bump'). As mentioned in Ref. [24], this is not 
uncommon since mechanisms that produce a sharp change in amplitude (such as required for the infrared 
cutoff) can also have an overshoot followed by damped oscillations [35]. Here, we focus on inflationary 
models with an infrared cutoff followed by a localized excess power (referred to as 'bump') in the power 
spectrum as a feasible possibility in explaining this deviation. This seems to be a promising avenue (such 
as [36]). Such features have also been seen by large-scale structure surveys [13, 34, 37, 38]. We have carried 
out an exhaustive parameter estimation with such a class of primordial spectra using the Bayesian approach 
by applying the Markov Chain Monte Carlo (MCMC) technique [39]. 

In the sections to follow we first describe briefly in section II, the various forms of primordial power 
spectra considered in this study. This is followed in section III by a brief overview of the methodology 
adopted (MCMC). In Section IV we give the parameter estimates with such a class of primordial power 
spectra which is followed by the final section on discussion and conclusion. 



While reconstructing the initial conditions of the universe it is important to devise a method that can 
efficiently constrain the shape of the power spectrum of primordial matter density fluctuations. It has 
been shown in [24] that the direct recovery of the primordial power spectrum from the CMB (WMAP) 
observations indicate a sharp infrared cutoff around the horizon scale (fc» ~ kh = 2n/rjo 7 where 770 is the 
present horizon) followed by a 'bump' [40]. The resulting angular spectrum gives a significantly better 
likelihood than a simple (scale-invariant) Harrison-Zel'dovich, scale-free power spectrum and also a simple 
monotonic infrared cutoff. This motivates us to consider a broader variety of primordial power spectra with 
infrared cutoffs at low wavenumber arising from physics in the initial phase of a (limited) inflationary epoch. 

In this section we describe such primordial power spectra with features which have been considered in this 
work and each model is motivated by some plausible inflationary scenario. Fig. 1 shows the best fit spectra 
for the class of primordial power spectra considered below. 



Before we go on to describe initial power spectra with features we would like to mention the commonly 
assumed power law power spectrum (referred to as 'Power Law' model here). This is generally parametrized 
in terms of the slow-roll approximation of the single-inflaton field and is given by 



Here n s (k) = [d\nV {k) / dink] + 1 is the scalar spectral index which is constant for power law models (n s = 1 
corresponds to a scale-invariant power spectrum) and ko is the scalar pivot scale. The power spectrum 
amplitude, A s determines the variance of the fluctuations. 

The various inflation motivated primordial power spectra models [41] with infrared cutoffs (with or without 
a 'bump') which have been explored in this work are described below. 



II. PRIMORDIAL POWER SPECTRA WITH FEATURES 



A. Scale-free Power Law model-PL(I) 




(2.1) 
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B. Exponential Cutoff model-EC(II) 



This model has a monotonic cutoff imposed on (multiplying) a scale-free power law spectrum. Such an 
infrared cutoff discussed and motivated by [30] leads to power spectrum of the form (referred to as 'Expo- 
cutoff' model here) 



V{k) = A s W 



1 



(2.2) 



In previous analyses using this model in [30, 31] the value of the parameter a was fixed (a = 3.35) to 
approximately mimic the sharpness of the cutoff in a specific scenario (See Kin.Dom model described later). 
In our work here we allow it to vary and obtain the best fit value. As discussed later, our results indicate 
that a larger value of a is preferred. 



C. Starobinksy model-SB(III) 



It was shown by Starobinsky [42] that if the effective scalar inflaton potential has a 'singularity' in the form 
of a sharp change in the slope (a 'kink'), it can create an infrared suppression in power spectrum of adiabatic 
perturbations at any chosen wavenumber fc* (referred to as 'Starobinsky' model here). The infrared cutoff 
is followed by the 'bump' that arises naturally as the first peak of a damped ringing. The effect of a kink 
(sharp change in the slope) at some point of the inflaton potential can be neatly expressed in terms of the 
analytic multiplicative transfer function applied on the underlying power spectrum V (k) as 



V(k) = V (k) V 2 (y,R*). (2.3) 

In the simpler form (SB) we restrict V {k) to be a simple power law. However, in general, V (k) can be of 
any form allowed by models of scalar field driven inflation. For example, in the later section we take it to 
be of an exponential form (Expo-Staro model). The transfer function 

V 2 (y,R*) = [l-3(R*-l)-((l--^)sm2y+-cos2y) + ^(.R„ - 1) 2 1(1 + 1) x 

y yZ y 2 ^ ^ 

(1 + l + (1 ^ )cos% _£ 8i „ %)] , (2 . 4) 

y 2 y 2 y 

with y = fc/fc». i?» is the ratio of the slope dV/d<j) which dictates the shape of the power spectrum, V{<j)) 
being the inflaton potential of the scalar field cf). The power spectrum V{k) in eq. 2.3 has a step-up like 
feature (going to larger k) for R r < 1 and a step-down like feature for i?» > 1. An infrared cutoff is created 
when R* < 1. This shape of the primordial power spectrum with a Starobinsky step not only mimics the 
sharp infrared cutoff but also subsequently produces the required localized excess in power (the 'bump'). We 
have carried out MCMC analysis and obtained best fit values for k* and 



D. Pre-inflationary Kinetic Domination model-KD(IV) 



The observable inflationary epoch could be preceded by a period of fast rolling of the inflaton field <p leading 
to pre-inflationary phase of kinetic domination. The difference of the vacuum in the kinetic domination (fast 
rolling) phase relative to the inflationary phase would imprint a feature in the power spectrum at large 
scales corresponding to first modes that crossed out of the Hubble radius at the onset of inflation [30] . The 
feature is an infrared cutoff akin to that first shown by Vilenkin and Ford [43] for a radiation dominated 
pre-inflationary phase. If the scale corresponding to the current horizon exited the horizon very soon after 
the onset of inflation, then the feature could explain the observed suppression of power at the horizon scale. 
The form of primordial perturbations is given by (referred to as 'Kin. Dom.' model here) 

V(k) = ^k\A-B\ 2 (2.5) 

where 
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FIG. 1: The class of primordial power spectra with features which have been used in this work are plotted. Here 
kh = 2-7T jr\ « 4.5 x 10~ 4 Mpc -1 , is the wavenumber corresponding to the Horizon scale for best fit ACDM model. For 
comparison we also give the power spectrum recovered from WMAP data by direct deconvolution [24] ('Recovered') 
and wavelet smoothed version that retains the most prominent features ('Recovered(sm)') (adapted from R. P. 
Manimaran et ai, in preparation). Notably, none of the model spectra are able to match up well to the form 
of infrared cutoff indicated by the direct deconvolution leaving room for further search of more likely inflationary 
scenario. 



and 



Jk/H 



B = 



\J 32 H J 7T 



2H 



H 



(2) 



k 

2H 



(2) 

Here H denotes the (physical) Hubble parameter during inflation while TLq and denote the Hankel 
function of the second kind with order and 1, respectively (for details refer to [30] and references therein). 



/(2) 



E. Pre-inflationary Radiation domination model- VF(V) 



For a pre-inflationary radiation dominated epoch the power spectrum was given by Vilenkin and Ford [43] 
as follows (referred to as 'VF' model here) 



P(k) = A s k 1 ^ — |e- 2 ^(l + 2iy) ~l-2y 2 \\ (2.6) 

where y — fc/fc*. The VF power spectrum can also provide an infrared cutoff with the required 'bump' after 
it. The infrared cutoff here is not as sharp as the one arising from kinetic domination in the pre-inflationary 
phase and is found to be relatively disfavored in our work. The cutoff scale fc* is set by the Hubble parameter 
at the onset of inflation. Here too, the current horizon scale crosses the Hubble radius very close to the onset 
of inflation. 
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F. Exponential-Starobinsky model-ES(VI) 



Encouraged by the reasonable success of the class of primordial spectra described in the earlier subsections 
we have undertaken a power estimation with yet another form of initial spectra tailored to mimic the 
cutoff and the 'bump' in the primordial power spectrum. As discussed earlier, the Starobinsky feature is a 
transfer function that modulates an underlying power spectrum expected from the inflaton potential without 
singularity and affects any inflationary scenario based on a inflaton potential with sharp change in slope 
(eq. 2.4). We imprint a Starobinsky break on an exponential cutoff spectrum (referred to as 'Expo-staro' 
model here) mainly motivated by the fact that these models provide a sharp reduction in large wavelength 
power needed to suppress the low multipolcs of the CMB. The Expo-staro model has both the exponential 
cutoff providing an adjustable sharpness of the cutoff as well the bump of the Starobinsky feature and is 
more fine tuned as reflected by the additional parameter, e*. This can be described as 



where V 2 {y,R st ) is the transfer function of the Starobinsky feature described by eq.( 2.4) and sets the 
ratio of the two cutoff scales involved (exponential cutoff and Starobinsky feature) . 

Using the recent measurement of the angular power spectrum of CMB anisotropy we undertake a full 
blown multi-parameter joint estimation of the parameters characterizing the infrared cutoff such as a 
and i?* (described above) together with the standard set of cosmological parameters. In the next section we 
briefly describe the parameter estimation method used in this work. 



Theoretical models for the determination of the power spectrum, which is responsible for structure forma- 
tion, starts from the assumption of a primordial V {k) at a sufficiently high redshift z 3> z eq , the rcdshift 
at matter-radiation equality. The primordial power spectrum is related to the linear CMB anisotropy power 
spectrum C ; Ml2 by the transfer function T, Ml2 (fc) such that [21] 



where Mi and Mi are the indices for the temperature and the two polarization modes of the power spectra. 

We use the well-known Monte Carlo Markov Chain (MCMC) technique [39, 44] to constrain the parameters 
with the observed data. Given a vector of cosmological data d (for instance, d can be the vector containing 
the 899 WMAP measurements of the temperature power spectrum Ci for I = 2 to 900) from which we want to 
measure a vector of cosmological parameters p = (r, f2 m , Q^, fc*,...). The cosmological parameterization 
has been carried out in terms of the following parameters, namely, the baryon density flbh 2 , cold dark matter 
density, Q c d m h 2 7 Hubble constant H a [45], and opacity to the surface of last scattering, r [46]. Apart from 
these, we include the parameters characterizing the inflationary primordial perturbations with infrared cutoff 
features as described in the previous section (fc», a, i?») along with the scalar spectral index of perturbations, 
n s [47] and the amplitude log[I0 10 A s ]. We consider a flat universe with Slk — 0. Table I gives the priors used 
in the study for the various inflationary and cosmological parameters. The ranges of these uniform priors 
have been varied to ensure that there is no dependence of priors on our results. We restrict our study to 
a cosmological flat model where the cosmological constant, 17a = 1 — Q c dm — fit, an d do not consider other 
forms of dark energy. The theory p is connected to the data d by the likelihood function, £(p, d), which 
gives the probability distribution for observing different d given a theoretical model p which in our case is 
the ACDM model. In Bayesian analysis, one inserts the actual observed data and and reinterprets C(p, d) as 
an un-normalized probability distribution over the cosmological parameters p, optionally after multiplication 
by a probability distribution, which has the prior information. To place constraints on any single parameter, 
one marginalizes (integrates) over all the others. In the MCMC method a large set of points p,,i = I, 
. . . n, i.e. a Markov chain, is generated by a stochastic procedure such that the points have a probability 
distribution £(p, d). The basic MCMC algorithms used here are based on the Metropolis-Hastings and slice 
sampling methods. According to Metropolis-Hastings rule [49] higher likelihood trials are always accepted 
i.e. for £(p*, d) > £(p,, d), otherwise rejected and accepted only with probability £(p*, d)/£(pj, d) where 
p* is a new trial random point in the parameter space. 

The other method used in the implementation of MCMC is slice sampling [39]. In slice sampling [50] 
similar to the Metropolis method the target density p*(x) is evaluated at any point x. Here is a brief note 



V(k) = V (k) l-e-(*. V 2 (y,R«) 



(2.7) 



III. METHODOLOGY 




(3.1) 



G 



TABLE I: Uniform prior ranges for the parameters as used in the joint parameter estimation study. 

Parameter Prior (range) 

Baryon Density Q b h 2 0.005-0.1 

Dark Matter Density S! ira /i 2 0.01-0.99 

Angular size of Acoustic Horizon 9 0.5-10.0 

Optical Depth r 0.01-0.8 

Scalar Spectral Index n s 0.5-1.5 

Infrared Cutoff (wavenumber) k, (x lO^Mpc" 1 0.00001-0.0009 

Expo Cutoff Steepness Parameter a 0.1-10.0 

Starobinsky Parameter R, 0.01-0.99 

Amplitude Parameter 6 logl0 10 A s 2.7-4.0 



"Acoustic horizon as defined in [45] 

b Amplitude parameter as defined in [48] 



on slice sampling algorithm. Consider a one-dimensional slice sampling method for making transitions from 
a two-dimensional point (x,u) lying under the curve p(x) to another point (x',u') lying under the same 
curve, such that the probability distribution of (a;, u) tends to a uniform distribution over the area under the 
curve p*(x), whatever initial point we start from. A single transition (x,u) — > (x',u') of a one-dimensional 
slice sampling algorithm has the following basic steps. A vertical coordinate u' and an horizontal interval 
(x a ,Xb) enclosing x (i.e. a slice) is then drawn under the curve p*(x). Now x' is drawn between this interval 
and p*(a/) is evaluated. If p*(a/) > u' then accept it otherwise choose a different interval and continue. 
There can be slight variations to this, for instance, 'stepping out' or 'shrinking' methods can be used for 
choosing the interval (x a ,Xb), the distribution under the curve always remaining uniform. Slice sampling is 
more robust to the choice of parameters like step sizes which is self tuning in this method and hence, we get 
better estimates when implementing MCMC with the new parameter set. However, slice sampling is slightly 
slower in terms of computer run-time. 

We have evaluated the likelihood function C and obtained the best fit from the parameter searches using 
a modified version of CosmoMC [39]. We typically ran multiple chains on a 64 bit 8x4 Alpha-Server ES45 
68/1250 System (1.25 GHz each) to optimize our parameter set. We have first conducted the analysis 
with only WMAP TT data wherein the discrepancy in power is more apparent. Later on the analysis 
has been extended to include both the WMAP TT and TE angular power spectra along with the CBI 
(Cosmic Background Imager), ACBAR (Arcminute Cosmology Bolometer Array Receiver) and VSA (Very 
Small Array) datasets of the CMB for small scales (I > 700). These external datasets augment the quality 
of the analysis as well as help in alleviating parameter degeneracies. For instance, these are essential to 
break parameter degeneracy between the inflationary parameter, n s and optical depth, r. In addition, high 
resolution CMB datasets also help to tightly constrain the parameters for example, fl m , amplitude A s , the 
running in the spectral index (not used here) [51-53]. 



IV. PARAMETER ESTIMATES 



We have carried out a complete Bayesian analysis by varying and marginalizing over all the parameters 
considered in this work to obtain the best fit values. Table II and Table III sums the results of the MCMC 
simulations with the WMAP data only and WMAP data plus the addition of the external datasets, respec- 
tively Tables II and III mention the values of some of the relevant parameters (the search has been made for a 
larger parameter space) along with their l-cr error bars and their goodness-of-fit parameter (Xe// = — 2 ln£) 
values for the different models. 

The best likelihood value (with the WMAP TT data) is given by ESa-b(VI) whereas EC(II) and SB (III) 
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TABLE II: Best fit values of parameters specifying the initial power spectrum (fc„, a, R*, n s ) and other relevant 
cosmological parameters for a class of model power spectra with a infrared cutoff (dataset used: WMAP TT data). 
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te» = 0.2 
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"Rcionization epoch as defined in [46] 

likelihood values are comparable. We obtain e* = 0.2 (ES-a(VI)) by carrying out MCMC with ES(VI) model 
which gives the best likelihood for this model while e* = 0.75 (ES-b (VI)) gives the best shape (i.e. closest 
to the recovered power spectrum) but is slightly less likely. In Fig. 2 we have plotted the angular power 
spectrum C; obtained using the best fit values after the MCMC analysis. We corroborate the preferred 
value for fc» ~ 3 x 10~ 4 Mpc -1 as obtained by [31] for the exponential cutoff power spectrum (EC(II)). An 
important difference to be noted is in the value of a. We obtain better likelihoods as the value of a increases. 
In fact, it can be seen from Table II that the best likelihood is obtained for a steeper exponential cutoff, i.e. 
a = 9.6 (— ln£ = 485.2 for a = 3.35) in contrast to earlier studies [31]. This is also in good agreement with 
the recovered power spectrum of [24] where again a very sharp cutoff (a = 10), gave a better fit to WMAP 
data. 

Fig. 3 shows the mean likelihoods for fc* for the models using WMAP TT data. We note that compared 
to other models the probability for model ES-a(VI) is sharply peaked on fc* ~ 3.2 x 10~ 4 Mpc -1 . This 
however, does not translate to tight \-a error bars in the corresponding table due to the bimodal probability 
distribution with a smaller, broad secondary peak at low wavenumber for this model. 

The values of the spectral index, n s and cosmological constant, Q\ affect the shape and the amplitude of 
the power spectrum of matter distribution. We would like to mention that getting the estimate of r from 
the low multipoles is important to break parameter degeneracies. As more parameters are added, precision 
continues to reduce (unless new degeneracies are broken) as many of the parameters are correlated. Deviation 
in primordial spectral index n s is a consequence of strong correlation with r [54]. In Fig. 4 we show the 
correlation between the optical depth r and the inflationary parameters, a, i?* along with the correlated 
n s . We notice that the increase in r and fc* is compensated by a corresponding increase in n s . Therefore, 
we find that inflationary models such as SB (III), KD(IV), and ESa-b(VI) models (with cutoff plus 'bump') 
provide better fit to the observed WMAP data. Another interesting result that emerges is the difference 
in r inferred with and without the WMAP polarization TE data. It is much less marked for models with 
infrared cutoff than for the pure Power Law model. Perhaps this could be an indirect evidence in favour 
of an infrared cutoff. We consider the VF model as not a very good fit to the data and leave it out of 
any further discussion. Table III gives the parameter search results with the external CMB datasets (CBI, 
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FIG. 2: The angular power spectra (C;) for the different primordial power models along with their best fit likelihood 
values are given in the figure. The plot also shows the WMAP binned C;S for comparison 




FIG. 3: The plot shows the mean likelihoods of the cutoff parameter, fc* obtained after performing MCMC for the 
models as discussed in the text using WMAP TT data. Notably, the probability for model ES-a(VI) is sharply peaked 
around fc» ~ 3.2 x 10" 4 Mpc" 1 . 
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TABLE III: Best fit values of parameters specifying the initial power spectrum (k*,a, R*,n s ) and other relevant 
cosmological parameters obtained using a class of model power spectra with a infrared cutoff (dataset used: WMAP 
(TT+TE) + External (CBI, ACBAR, VSA) CMB datasets) 



Parameter Expo-cutoff Starobinsky Kin. Dom. VF Expo-staro(a) ' Expo-staro(b)* Power Law 
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ACBAR, VSA). In this case too we report the best likelihood with ESa-b(VI) P(k) compared to the rest 
of the pack. To investigate the effect of priors we performed a number of MCMC simulations with different 
priors within the ranges given in Table I and found that the dependence of the results on priors is weak [55]. 

We note that introducing a cutoff fc* and a 'bump' at low wavenumber in the primordial power spectrum 
provides a better fit to the observed CMB data. Such scenarios can explain the lack of power at low multipoles 
as shown in Fig. 2. 

V. DISCUSSION AND CONCLUSION 

The cosmological model and parameter constraints deduced from observations are sensitive to assump- 
tions regarding the power spectrum of primordial density perturbations. The suppression of power in the 
quadrupole of the CMB anisotropy could be signaling a radical departure (in the form of an infrared cutoff) 
from the near-scale invariance of the primordial spectrum predicted by simplest inflationary scenarios. The 
very presence of a feature in the primordial spectrum indicates modification to the simplest inflationary 
scenarios and calls for an exploration of extensions based on plausible physical mechanisms. 

We show that sharp features, such as a infrared cutoff near the horizon scale and a 'bump', in the primordial 
power spectrum provides a likely explanation to power suppression at low multipoles as seen in the recent 
WMAP data. Although a steep monotonic cutoff tends to pull down the power in the next few higher 
multipoles above the quadrupole and octopole but this is well compensated by the sudden excess power 
following the sharp rise. This allows a steep cutoff to match the low quadrupole and octopole without 
suppressing the higher multipoles. This is evident in the studies where the primordial power spectra was 
directly deconvolved from the data [24, 34] . In the light of the analysis conducted in this work we give more 
credence to inflationary models with cutoffs than put forth in earlier studies [30, 31]. It can be seen through 
this study that a break in the initial power spectrum confronts the CMB data better. This is also borne 
out in recent Bayesian evidence assessment of a different set of model spectra with infrared cutoff [29] . The 
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FIG. 4: The plots show the posterior distribution of the parameter r plotted against the parameters k* and a for 
exponential cutoff EC(II) model (Top left panel), fc* and R, for Starobinksy SB(III) model (Top right panel) and k*, 
a and R* for expo-staro ES-a(VI) model (bottom panel). Points are colored according to the value of n 3 . 

question however, remains of the location of this cutoff. Our study suggests that a cutoff of fe* ~ 3 x 10~ 4 
Mpc -1 with WMAP TT data and 1.1 x 10~ 4 < K < 9.0 x 10~ 4 Mpc -1 (inclusive of error bars) with all 
the CMB datasets. However, more detailed and wider exploration of the infrared cutoff mechanism within 
inflation is required. A comparison of the best fit model spectra with the recovered spectra shows that the 
form of infrared cutoff has substantial scope for improvement. This keeps open the pursuit of (physically 
motivated) mechanisms within simple inflation that produces the form of infrared cutoff suggested by the 
data. 

Studies such as [36] also assert that the presence of such features if confirmed further by observations 
would be attributed to the initial power spectrum and not due to oscillations produced by any other means 
as suggested in Ref. [56] . Other feasible possibilities which produce such features in the initial power spectrum 
are the multiple-field inflationary models. Some of these have been tested with the recent observed CMB 
data such as the double inflation [33] scenario. However, such studies indicate that the present CMB data 
(WMAP first year data) strongly favour single-inflaton field models. Further discussion on the classes of 
inflationary models supported by observational data can only be made and asserted after future advancement 
in measurements, in particular, those of CMB polarization. The polarization power spectra have different 
parameter dependence as opposed to the temperature power spectra. This aids in breaking some parameter 
degeneracies. Moreover, since their source of origin is different and the polarization spectra directly probes 
the epoch of last scattering so polarization can be a better probe of the inflationary epoch. 
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